Constrained Opinion Dynamics: Freezing and Slow Evolution 
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We study opinion formation in a population of leftists, centrists, and rightist. In an interaction 
between neighboring agents, a centrist and a leftist can become both centrists or leftists (and 
similarly for a centrist and a rightist), while leftists and rightists do not affect each other. For 
any spatial dimension the final state is either consensus (of one of three possible opinions), or a 
frozen population of leftists and rightists. In one dimension, the opinion evolution is mapped onto 
a constrained spin-1 Ising model with zero-temperature Glauber kinetics. The approach to the 
final state is governed by a t~^ long-time tail, with ip a non-universal exponent that depends on 
the initial densities. In the frozen state, the length distribution of single-opinion domains has an 
algebraic small-size tail a;-^''-'*) with average domain length L^^ , where L is the length of the 
system. 

PACS numbers: 64.60.My, 05.40.-a, 05.50.-|-q, 75.40.Gb 



One of the basic issues in opinion dynamics is to un- 
derstand the conditions under which consensus or diver- 
sity is reached from an initial population of individuals 
(agents) with different opinions. Models for such evolu- 
tion are typically based on each agent freely adopting a 
new state in response to opinions in a local neighborhood 
[01 . The attribute of incompatibility - in which agents 
with sufficiently disparate opinions do not interact - has 
recently been found to prevent ultimate consensus from 
being reached H, |j. Related phenomenology arises in 
the Axelrod model ||, a simple model for the forma- 
tion and evolution of cultural domains. The goal of the 
present paper is to investigate the role of incompatibility 
within a minimal model for opinion dynamics. This con- 
straint has a profound effect on the nature of the final 
state. Moreover, there is anomalously slow relaxation 
to the final state. While we primarily frame our dis- 
cussion in terms of opinion dynamics, our results apply 
equally well to the coarsening of spin systems. In the 
latter context, we obtain a new non- universal kinetic ex- 
ponent in one dimension that originates from topological 
constraints on the arrangement of spins. 

We consider a ternary system in which each agent can 
adopt the opinions of leftist, centrist, and rightist. The 
agents populate a lattice and in a single microscopic event 
an agent adopts the opinion of a randomly-chosen neigh- 
bor, but with the crucial proviso that that leftists and 
rightists are considered to be so incompatible that they 
do not interact. While a leftist cannot directly become 
a rightist (and vice versa) in a single step, the indirect 
evolution leftist =J> centrist =J> rightist is possible. Our 
model is similar to the classical voter model |^ and also 
turns out to be isomorphic to the 2-trait 2-state Axelrod 
model P, 0| . Due to the incompatibility constraint in our 



model, the final opinion outcome can be either consensus 
or a frozen mixture of extremists with no centrists. Fig- 
ure |l| shows a typical frozen state on the square lattice 
(with periodic boundary conditions). Notice the nested 
enclaves of opposite opinions and the clearly visible clus- 
tering. 
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FIG. 1: Typical frozen final state in our opinion dynamics 
model on a 100 x 100 square lattice for po = 0.1. The two 
extreme opinions are represented by black and white squares. 

We can exploit the connection between the voter model 
and our opinion dynamics model to infer the final state 
of the latter. If we temporarily disregard the difference 
between leftists and rightists, the resulting binary system 
of centrists and extremists reduces to the voter model, for 
which one of two absorbing states — either all centrists 
or all extremists — is eventually reached. In the context 
of the ternary opinion system, the latter event can mean 
either a consensus of extremists or a frozen mixed state 
of leftists and rightists, as depicted in Fig. |l|. 

Because of the underlying voter model dynamics, the 
average density of each species is globally conserved in 
any spatial dimension. Therefore {pi(i)) — pi{t — 0), 
where i refers to one of the states (-I-, 0, — ) and the an- 
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gle brackets denote an average over all realizations of the 
dynamics and over all initial states that are compatible 
with the specified densities. As a result of this conserva- 
tion law, with probability Pq = po the final state consists 
of all centrists and with probability 1 — po there are no 
centrists in the final state. In the latter case, there can 
be either a consensus of + (this occurs with probability 
P+), consensus of — (probability P_), or a frozen mixed 

state (probability ). Figure ^ shows the dependence 

of these final state probabilities on po in the mean-field 
limit (where all agents are interconnected) in the sym- 
metric case (0+ = (0- = (1 — /3o)/2. In one and two di- 
mensions, the final state probabilities are nearly identical 
to the mean- field predictions when p+ = but differ- 
ences become apparent in the strongly asymmetric cases 
of p+ » P-. The final state probabilities also depend 
very weakly on the system size. 
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FIG. 2: Probability for the occurrence of a given final state 
as a function of po for p+ = p-. Here P+ is the probability 
that + consensus is reached and P+_ is the probability that 
the final state is a frozen mixture of -I- and — . 

We now focus on the one dimensional case. Here our 
opinion dynamics model is equivalent to a constrained 
spin-1 Ising chain that is endowed with single-spin fiip 
zero-temperature Glauber kinetics 0, with leftist, cen- 
trist, and rightist opinions equivalent to the spin states 
— , 0, and respectively. The incompatibility constraint 
means that neighboring -I- and — spins do not interact. 

This Ising model picture suggests that the best way to 
analyze the dynamics in one dimension is to reformulate 
the system in terms of domain walls. There are three 
types of domain walls: freely diffusing mobile domain 
walls between -1-0 and between —0, denoted by M+ and 
M_ , respectively, and stationary domain walls 5* between 
H — . The mobile walls evolve by 



M± + M± 



A/± + S. 



(1) 



M±+S^ A/zr 



(2) 



Thus stationary domain walls are dynamically invisible; 
their only effect is that the sign of a mobile wall changes 
whenever it meets a stationary wall, after which the latter 
disappears (Fig. |^). The inertness of the stationary walls 
is reminiscent of kinetic constraints in models of glassy 
relaxation. These constraints typically lead to extremely 
slow kinetics |, | 0, as is also observed in our 
opinion dynamics model. 

-+-+-0+ -+ +- 




FIG. 3: Space-time representation of the domain wall dynam- 
ics. Time runs vertically downward. The spin state of the 
domains and the identity of each domain wall are indicated. 

The rate equations corresponding to the processes in 
Eqs. (1) and ^) are 



S = + 



(3) 



When a mobile wall hits a stationary wall, the former 
changes its sign while the latter is eliminated via the 



These give the asymptotic behaviors M (x and S oc 
^-1/2^ Thus an approximate rate equation approach al- 
ready predicts that stationary domain walls decay slower 
than mobile walls. 

An important subtlety in the arrangement of domain 
walls is that an arbitrary initial opinion state necessarily 
leads to an even number of mobile walls between each 
pair of stationary walls. It is also easy to verify that do- 
main wall sequences of the form . . . . . . cannot 

arise from an underlying opinion state. These topological 
constraints play a crucial role in the kinetics. 

The exact density of mobile walls can be obtained by 
again mapping the constrained spin-1 system onto a spin- 
1/2 system that is equivalent to the voter model. In this 
mapping, we consider both -I- and — spins as compris- 
ing the same (non-zero) spin state, while the zero spins 
comprise the other state. With this identification, the 
reduced model is just the spin-1/2 ferromagnetic Ising 
chain with zero-temperature Glauber dynamics and no 
kinetic constraint. In this reduced model, domain walls 
M+ and Af_ are indistinguishable and they diffuse and 
annihilate when upon colliding. The density of mobile 
walls M(t) = M+{t) + M-{t) is known exactly for arbi- 
trary initial conditions from the original Glauber solution 
[Q. For initially uncorrelated opinions and if the mag- 
netization of the spin-1/2 system - here the difference 
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between the density of non-zero and zero spins - equals 
mo, then |p| 



M{t) = 



1 — m, 



°e-2* [/o(2t)+/i(2i)] 



(4) 



where Ik is the modified Bessel function of index k. 

In the spin-1 system, mg = /9+ + p_ — po, or mg = 
1 — 2pQ by normalization. If the initial densities of + and 
— opinions are equal, then M+(i) = M_(i) and their 
densities are 

M±{t) = pail - po) e-^' [Io{2t) + h{2t)] 

^ po(l-po)(^t)-^/'. (5) 

As expected, the mobile wall density asymptotically de- 
cays as because of the underlying diffusive dynam- 
ics. However, we find numerically that the density of 
stationary domain walls S{t) decays as 



(6) 



with a non-universal exponent tpiPo) that goes to zero 
as po ^ (Fig. I). 
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FIG. 4: Stationary domain walls density versus time on a 
double logarithmic scale for the initial conditions po = 0.02, 
0.04, 0.10, 0.20, and 0.40 (top to bottom). The respective ex- 
ponent estimates are 0.013, 0.026, 0.065, 0.13, and 0.29. Data 
are based on 100 realizations on a 5 x 10^-site chain. Inset: 
Stationary domain wall density for initially uncorrelated walls 
for Po = 0.02, 0.10, and 0.40 (top to bottom). The solid line 
has slope —3/8. 

To help understand the mechanism for the slow decay 
of the stationary domain wall density, we also simulated 
a test system with spatially uncorrelated domain walls. 
While such a domain wall state cannot arise from any 
initial set of opinions, we can prepare directly an un- 
correlated arrangement of domain walls with prescribed 



densities. For any initial condition in this test system, the 
stationary wall density decays as t~^/^ (inset to Fig. H), 
consistent with known results on persistence ^ . 

Here persistence refers to the probability that a given 
lattice site is not hit by any diffusing domain wall. For 
the kinetic spin-1/2 Ising model, the persistence proba- 
bility decays as t^^ , with 6 — 3/8 ||T^, independent of the 
initial domain wall density, when the walls are initially 
uncorrelated. Thus the topological constraints imposed 
on the domain wall arrangement by the initial opinion 
state appear to control the dynamics. 

These topological constraints lead to the initial- 
condition dependence of the amplitude in the mobile wall 
density (Eq. (|)). This arises because for po ^ the sys- 
tem initially consists of long strings of stationary walls 
that are interspersed by pairs of mobile walls, and their 
survival probability is proportional to their initial (unit) 
separation , leading to the asymptotic density for mo- 
bile walls is M ~ 2po/V7rt (Eq. (g)). We now exploit this 
observation to estimate the density of stationary walls as 
Po — > 0. Within a rate-equation approximation, the den- 
sity of stationary domain walls decays according to 



S^-kMS. 



(7) 



While such an equation is generally inapplicable in low 
spatial dimension, we can adapt it to one dimension 
by employing an effective time-dependent reaction rate 
k ^ ^Ijixt |T^. This is just the time-dependent 
fiux to an absorbing point due to a uniform initial back- 
ground of diffusing particles; such a rate phenomenolog- 
ically accounts for effects of spatial fluctuations in one 
dimension. Substituting the asymptotic expression for 
M (t) from Eq. (^) and the reaction rate k ~ y/2/TTt into 
this rate equation, we find that the density of stationary 
walls decays as t'"^ with tpipo) = VSpq/tt as po — > 0. It 
is the amplitude in the density of mobile domain walls 
that ultimately causes the slow decay in the stationary 
wall density. 

A more compelling way to determine ipiPo) is via per- 
sistence in the g-state Potts model. Because the ini- 
tial magnetization in the reduced spin-1/2 system is 
niQ = 1 — 2po, it has been argued (see e.g. ||l8[) that 
this system should be identified with the g-state Potts 
model with mo — 2/q — 1, or g = (1 — po)~^- Using the 
exact persistence exponent for the g-state Potts model 
with Glauber kinetics | [T^ , and identifying ip with this 
persistence exponent, we obtain 



8 TT^ 



cos 



l-2po 
V2 



(8) 



with the limiting behavior '>p{po) — > 2po/7r as po — > 0. 
This asymptotics agrees with our numerical results for 



Pa 



< 



0.4 (Fig. H) but then deviates for larger po, where 
ip{po) must monotonically approach 1/2 as po 1- It 
should be noted that in this identification with persis- 
tence, we have ignored the creation of stationary inter- 
faces due to the meeting of mobile domain walls. This 
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FIG. 5: Comparison of the exponent i/) from Eq. and from 
the simulation data of Fig. ^ (circles). 

An important characteristic of the system is the spatial 
arrangement of domain walls. From our simulations, the 
mean distances between nearest-neighbor MM and MS 
walls both appear to grow as t^/'^ due to the diffusive mo- 
tion of mobile domain walls. The distributions of these 
two distances both obey scaling, with scaling function of 
the form ze~^ , where z = x{t)/{x[t)) is the scaled sepa- 
ration between walls. In contrast, the distances between 
neighboring stationary walls xgs appear to be character- 
ized by two length scales. There are large gaps of length 
of the order of t^/'^ that are cleared out by mobile walls 
before they annihilate, but there are also many very short 
distances remaining from the initial state (Fig. |3|). The 
corresponding moments {xgg{t))^/'^ reflect this multiplic- 
ity of scales, with {xgg{t))^^'' approaching a t^/^ growth 
law as fc — > oo, and growing extremely slowly in time for 

Finally, we quantify the frozen final state by the mag- 
netization distribution P{m), namely, the density differ- 
ence between -I- and — spins. This distribution becomes 
broader as po increases (Fig. reflecting the fact that 
there is progressively more evolution before the system 
ultimately freezes. For small po, P{m) has a tail. 
We may explain this result by considering the evolution of 
a single pair of mobile walls. This pair annihilates at time 
t with probability density n(i) cx f'^l'^ . The total mag- 
netization of the resulting frozen state scales as t^l"^ since 
the domain wall pair annihilates at a distance x cx t^/^ 
from its starting point. Then from P{m) dm — n(i) dt, 
together with n(i) cx t^'^/^ and to cx pot^^^, we obtain 
P{m) cx pQin^^. While the argument has been formu- 
lated in one dimension, we expect this result to apply in 
all spatial dimension. 

The frozen state is reached when t — Tf oc L^; this 
is the time needed for mobile domain walls to diffuse 
throughout the system and thus be eliminated. At this 
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FIG. 6: Magnetization distribution P{m), (m — p+ ~ p-) in 
the frozen final state for po ~ 0.02 (10^ realizations), 0.04, 0.1, 
and 0.2 (10* realizations) on a 5000 site linear chain. Inset: 
P(m) for Po = 0.02 on a double logarithmic scale to illustrate 
the tail. The straight line has slope —2. 
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FIG. 7: Domain length distribution F{x) in the frozen final 
state for the same data in Fig. ^. The data have been binned 
over a small logarithmic range to reduce fluctuations. Tabu- 
lated are the slopes from each data set and the expected value 
2(1 — ■(/)) (parentheses) from our simulation result for t/j. 



time, the density of stationary domain walls is of the 
order of S" oc T^'^ cx L~^'^. Thus the average length of 
single-opinion domains is (x) cx L'^'^ . Numerically we find 
that the domain length distribution has a power law tail, 
F{x) cx x^^, with 1 < /Li < 2. The lower bound ensures 
normalizability while the upper bound implies that (x) 
diverges as i — > cx3. From the above power law form, 
(x) = J dxxF{x) cx i^^^. This matches our previous 
estimate of {x) ^ L^^ when /i = 2{1 ~ ip). Figure |^ 
shows the length distribution of single-opinion domains 
in the frozen state. Direct estimates of the exponent fi 
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from this plot are in good agreement with the exponent 
relation /i = 2(1 — ip), with tp obtained from the time 
dependence of the mobile wall density in Fig. ^. 

In summary, the constraint that extremists with op- 
posite opinions cannot influence each other substantially 
slows down opinion dynamics. In one dimension, the den- 
sity of stationary interfaces between neighboring -I- and 
— spins decays as t~'^ , with ip{po) ~ "^Po/tt as po — > 0. 
This slow dynamics is a consequence of the subtle topo- 
logical constraints on the arrangement of domain walls, 
together with the incompatibility constraint that neigh- 
boring -|- and — spins do not interact. 

The final opinion outcome depends non-trivially on the 
initial densities of the leftist, rightist and centrist opinion 
states. With probability po the final population consists 



of only centrists, while with probability 1 — po the final 
population does not contain any centrists. Analytically, 
we only know that consensus of centrists is reached with 
probability Pq — po; the determination of the comple- 
mentary final state probabilities for general initial condi- 
tions P-f, P_, and P-^ remains an open challenge. Fi- 
nally, a frozen mixture final state would not exist if there 
is a non-zero rate (however small) for opposite extrem- 
ists to influence each other directly. Perhaps this lack of 
direct influence can account for the sad phenomenon of 
the proximity of incompatible populations in too many 
locations around the world. 

We thank S. Majumdar for helpful correspondence. We 
are also grateful to NSF grant DMR9978902 for partial 
support of this work. 
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